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I.  INTRODUCTION 

The  geometrical  theory  of  diffraction  (GTD)  [Keller, 1953, 1956, 1962] 
provides  a very  convenient  and  accurate  method  for  analyzing  the  problem 
of  high  frequency  plane  wave  scattering  by  smooth  convex  surfaces  except 
in  the  vicinity  of  the  shadow  boundaries  and  caustics.  Thus,  approximate 
asymptotic  high  frequency  results  which  remain  valid  even  within  the 
regions  where  the  GTD  fails  are  developed  in  this  paper  for  the  case  of 
a smooth,  convex  cylinder  illuminated  by  a plane  wave.  This  work  is 
motivated  by  the  need  to  analyze  the  radiation  patterns  of  antennas 
mounted  near  curved  surfaces  such  as  an  aircraft  fuselage  [Burnside, 1972, 
Burnside,  et  al . ,1975]  or  a ship  mast. 

In  the  GTD  analysis  of  the  scattering  of  waves  by  an  impenetrable, 
smooth  convex  cylinder  [Keller, 1956],  the  total  exterior  field  is 
associated  with  the  usual  incident  and  reflected  rays  of  geometrical  optics, 
together  with  the  surface  diffracted  rays  introduced  by  Keller.  The 
geometric  optical  rays  do  not  penetrate  the  convex  body,  and  hence  their 
field  is  zero  within  the  shadow  region.  The  surface  diffracted  rays  are 
launched  by  the  incident  ray  which  grazes  the  surface  as  shown  in  Figure 
1;  these  surface  rays  propagate  into  the  shadow  region  and  entirely  account 
for  the  field  therein.  In  that  pure  ray  optical  analyses  fail  at  and  near 
shadow  boundaries,  the  GTD  solution  is  valid  in  regions  I and  III  as  shown 
in  Figure  1.  It  fails  in  the  penumbral  or  transition  region  adjacent  to 
the  shadow  boundary  (SB)  which  is  indicated  as  region  II  in  Figure  1. 

The  angular  extent  of  this  transition  region  is  of  the  order  (2/kpg(Q-|  ))'/3 
where  k refers  to  the  wavenumber  of  the  surrounding  medium  which  is 
assumed  here  to  be  free  space,  and  Pg(Ql ) is  the  radius  of  curvature 
of  the  cylinder  at  the  point  of  grazing  incidence,  i.e.,  at  Q-|.  GTD 
also  fails  in  regions  IV,  V,  and  VI  shown  in  Figure  1 which  are 
usually  referred  to  as  the  surface  or  caustic  boundary  layer  regions. 

Regions  IV  and  V are  close  to  the  surface  of  the  obstacle  which  is  a 
caustic  of  the  surface  diffracted  rays;  whereas,  region  VI  is  in  the 
neighborhood  of  Q-|  which  is  a caustic  of  the  reflected  ray  for  grazing 
incidence.  An  asymptotic  evaluation  of  the  field  within  the  shadow 
boundary  transition  regions,  the  surface  boundary  layer  regions  and  also 
their  common  regions  of  overlap,  is  complicated  by  the  fact  that  the 
field  must  change  rapidly  but  smoothly  from  one  form  to  another  across 
these  regions. 

The  problem  of  estimating  the  fields  within  the  transition  regions 
associated  with  the  diffraction  of  waves  by  smooth,  convex  surfaces 
has  received  much  attention,  especially  in  connection  with  the  theory 
of  radio  wave  propagation  around  the  earth.  Fock  made  significant 
contributions  to  this  area  by  developing  a general  asymptotic  theory 
for  the  diffraction  of  waves  by  large  convex  surfaces  [Fock, 1957].  Fock 
expressed  his  general  solution  to  this  problem  in  terms  of  a canonical 
integral.  In  particular,  Fock  considered  the  Fresnel  diffraction  by  a 
sphere  [Fock, 1951]  in  which  he  approximated  the  canonical  integral 
asymptotically  within  the  transition  region;  as  a result,  he  was  able 
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Figure  1.  Rays  associated  with  the  plane  wave  scattering 
by  a smooth  convex  cylinder. 

to  obtain  a simpler  solution  in  terms  of  tabulated,  universal  functions 
of  a single  parameter.  That  result  [Fock,1951]  is  valid  for  heights  of 
the  source  and  observation  points  above  the  sphere  which  are  small 
compared  to  the  sphere  radius.  In  another  important  paper  on  this 
subject  by  Wait  and  Conda  [1959],  this  condition  was  relaxed  so  that 
the  source  and  the  observer  could  now  be  far  from  the  diffracting 
surface.  The  results  in  [Wait  and  Conda, 1959]  are  in  terms  of  functions 
similar  to  those  in  [Fock,1951],  but  they  are  obtained  from  an  asymptotic 
analysis  of  the  canonical  problem  of  the  diffraction  of  a plane  wave  by 
a circular  cylinder.  However,  the  results  in  [Fock,1951]  and  [Wait  and 
Conda, 1959]  do  not  reduce  uniformly  to  the  usual  6TD  ray  solutions  for 
the  lit  and  shadow  regions  far  from  the  shadow  boundaries.  Consequently, 
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one  has  to  be  concerned  about  the  blending  of  the  transition  field  and 
the  ray  optical  field  solutions.  On  the  other  hand,  Logan  and  Yee  [1962] 
were  able  to  obtain  an  approximate  uniform  solution  which  recovers 
the  GTD  solution,  by  re-defining  the  parameters  in  Fock's  canonical 
integral  via  an  ad  hoc  procedure;  whereas,  Ivanov  [1971]  obtained  a 
uniform  result  also  in  terms  of  Fock's  canonical  integral  which  recovers 
the  GTD  solution,  by  employing  a coordinate  transformation  different 
from  that  of  Fock.  However,  both  of  those  uniform  solutions  are  far  too 
complex  for  numerical  computations  in  the  entire  transition  region.  It 
is  noted  that  the  results  in  [Fock, 1951]  and  [Wait  and  Conda,1959]  are 
valid  in  region  II  which  excludes  regions  IV  and  V of  Figure  1.  An 
extensive  bibliography  on  other  work  dealing  with  various  aspects  of 
this  diffraction  problem  may  be  found  in  [Logan  and  Yee, 1962;  Logan, 1959], 
and  [Borovikov  and  Kinber,1974], 

In  this  paper,  the  canonical  problem  of  plane  wave  scattering  by 
a perfectly  conducting  circular  cylinder  is  solved  in  a manner  analogous 
to  that  in  [Fock,1951j,  and  some  asymptotic  approximations  are  introduced 
heuristically  to  obtain  a uniform  result  for  the  total  field  which  is  valid 
within  the  transition  region  outside  the  surface  boundary  layer,  and  which 
automatically  recovers  the  usual  GTD  or  ray  optical  field  solution  exterior 
to  this  transition  region.  Hence,  the  present  solution  need  not  be  in- 
troduced separately  within  the  transition  region,  as  a correction  to  the 
GTD  solution.  Furthermore,  the  present  asymptotic  result  is  expressed  in 
the  convenient  format  of  the  GTD;  and  it  allows  a simple  interpretation 
in  terms  of  the  GTD  rays.  Only  the  leading  terms  in  the  asymptotic  ap- 
proximation are  retained  in  this  analysis.  This  uniform  result  is  expressed 
in  terms  of  tabulated,  universal  functions  which  are  similar  to  those  oc- 
curring in  [Fock, 1951]  and  [Wait  and  Conda,1959];  hence,  the  present  result 
is  very  convenient  for  engineering  applications.  The  analytical  details 
pertaining  to  the  solution  of  the  canonical  problem  are  presented  in 
section  II;  these  results  are  generalized  to  the  convex  cylinder  case  in 
section  III.  Since  the  uniform  result  presented  here  is  not  valid  in  the 
close  vicinity  of  the  surface  and  the  surface  boundary  layer,  a separate, 
approximate  representation  for  the  field  which  is  valid  very  close  to  and 
on  the  surface  is  obtained  in  section  IV.  Numerical  results  are  presented  in 
section  V to  indicate  the  accuracy  of  the  asymptotic  results  developed  in 
sections  III  and  IV. 


II.  DEVELOPMENT  OF  A UNIFORM  RESULT  FOR  THE 

CIRCULAR  CYLINDER  CASE 

An  asymptotic  high  frequency  analysis  of  the  canonical  problem  of 
electromagnetic  plane  wave  scattering  by  a perfectly  conducting  cir- 
cular cylinder  in  free  space  is  presented  in  this  section.  This 
analysis  leads  to  a uniform  asymptotic  result  for  the  field  in  the 
shadow  boundary  transition  region  which  does  not  include  the  close 
vicinity  of  the  cylinder.  The  geometry  of  this  problem  is  illustrated 
in  Figure  2.  Since  the  problem  is  two-dimensional,  it  reduces  to  a 
scalar  problem.  It  is  assumed  in  the  analysis  that  the  field  satisfies 
an  e-iw*  time  dependence  which  is  suppressed.  The  total  field  exterior 
to  the  cylinder  is  denoted  by  uT.  Thus,  u^  = u^  + us,  where  u1  denotes 
the  incident  field,  and  us  denotes  the  field  scattered  by  the  cylinder. 
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Figure  2.  Geometry  of  the  canonical  problem. 


If  the  incident  electric  field  has  only  a z component,  then  the 
problem  is  of  the  TMZ  (or  acoustic  soft)  type.  If  the  incident  magnetic 
field  has  only  a z component,  then  the  problem  is  of  the  TEZ  (or  acoustic 
hard)  type.  Let 

ui  = eJkx  = eJkPcos<!’  . (1) 


UT  satisfies  the  reduced  wave  equation  and  the  following  boundary  con- 
ditions. 


(v|  + k2)  uT  = 
Q uT(P,*)|p=a 


0;  = two  dimensional  Laplacian  operator. 


- 0 ; Q » 


1 for  TMZ  case, 
’jp-  for  TEZ  case. 


(2) 

(3) 
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us  satisfies  the  Sommerfeld  radiation  condition  as  p+“.  (4) 

It  is  noted  that  u1  and  us  satisfy  the  reduced  wave  equation.  A formal 
integral  solution  to  (2)  subject  to  (3)  and  (4)  is  given  in  [Bowman 
et.al.,1969]  as 

uT  = u + ucw  , (5a) 


where 

r"je  QJ  (ka)  ,9*  . , 

u = dv[Ju(kc)  - H<2)(kp)]  e-^ 

LU  v QHy2  (ka)  v 

or 


u 


j dvCH^dcp) 

_0°_  j £ 


QH^Cka) 

QH^(ka) 


Hj2)(kp)]e-Jv* 


(5b) 


(5c) 


with  e being  a positive  number,  however  small  and 
♦ = UI  - f-  » with  UI  15  * (or  U|  < j-  )• 


(5d) 


The  term  u in  (5a)  represents  the  higher  order  diffraction  effects 
which  become  negligible  for  large  ka;  more  will  be  said  about  this  term 
later  in  the  analysis.  Within  the  transition  region,  v * 0 (ka)  such  that 
one  may  employ  the  transformation 


= ka  + mx 


* Of  ■ 


(6a;6b) 


(1) 

4(2), 


and  replace  J (ka)  and  H'  '(ka)  in  (5b;5c)  in  the  transition  regions 
by 


(1) 

Jv( ka ) % (m/iT)“U(x);  H^(ka)  % +j(m/5T)  ^w^(t), 
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(7a;7b) 


where  the  Airy  functions  V(x)  and  w^x)  are  defined  as  [Logan  and  Yee,1962] 
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2j  V(x)  = w](x)-w2(x);  W-| ( x ) = Ip  [ dt 

O rv  Jy. 


.Tt-t3/3 
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The  approximation  in  (7)  is  valid  for  large  ka  and  v~0(ka).  It  is  assumed 

here  that  ka  is  large.  The  contour  r-|  runs  from  «e_j(2ir/3)  to  *-je,  and 
Tg  is  the  complex  conjugate  of  r-j.  in  (7),  terms  of  0(l/m2)  and  higher 
are  neglected.  Therefore,  within  the  transition  regions  ( 5b ; 5c ) become 
via  (7a;7b;7c  and  7d): 


u ^ m 


ld,[Jv(t) 


(kp ) 


♦ i fc) 


»-je 


Sfw9(x) 


H^)(kp)]e 


-j(ka+mx)^ 


(8a) 


or 


„ * m 
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l^H^)(kp) 


(0 


-*-je 


- H(?\(kp)]e-j(ka+mT)^ 

5k(x) 


(8b) 


with 


.11  , for  the  TM  case 

$=<  2 (8c) 

' -1  3 

I -m"  » for  the  TEZ  case 

The  field  u may  be  initially  decomposed  as  in  [Fock,1951]  for  the 
sphere  problem;  namely,  u is  expressed  as  the  sum  of  the  integral  in 
(8a)  over  the  contour  C2  from  0-je  to  ®-je,  and  the  integral  in  (8b) 
over  the  contour  c-|  from  -<»-je  to  0-je,  respectively.  After  rear- 
rangement of  terms,  and  the  use  of  H^Pjkp)  = 2J  , Jkp)-H^^(kp)  in 
the  integral  over  C2,  one  obtains,  V'T'  V^T'  V'T' 


o = I-|  + *2*  (9a) 


It  is  noted  that  I]  is  completely  independent  of  the  electrical  prop- 
erties of  the  cylinder.  One  may  approximate  the  first  integral  asymp- 
totically in  (9b)  via  the  method  of  steepest  descent;  it  is  seen  that  this 
term  constitutes  the  incident  field. 


r . , rJc  H<1)(kp)+H<2^(kp)  . 

j dT  Jv(t)(kp)e-Jv(t)*  - \ dv  a-  ; v e~J' 


cl+c2 


-*-je 


,-jkpsin.j/  = j kp  cos  1 4>  | ; 1 4>  | <k  (10) 


The  Hankel  functions  in  the  integrand  of  (10)  are  replaced  by  their 
Debye  approximations  prior  to  the  asymptotic  evaluation  of  the  integral. 
The  Debye  approximation  is  [Bowman  et.  al.,1969] 


a — J kp  SinY  + jvy  + T ; v=(cp  COSY 

0<ReY  <ir  . 


Since  0<ReY<ir  in  (II),  the  H^^(kp)  term  in  (10)  yields  u7  for  0<|cf>|<ir/2 
from  the  saddle  point  at  y^/^  - 1 <j> | , and  contributes  negligibly  else- 
where: likewise,  the  H(2)(kol  term  yields  u7  for  ir/2<UT<ir  from  the  sadc 


from  the  saddle  point  at  y=i r/2  and  contributes  negligibly  else- 

where; likewise,  the  H'2)(kp)  term  yields  u7  for  tt/2< I <t> | <tt  from  the  saddle 
point  at  y = | <j> | —ft/2 . vThe  remaining  integrals  in  (9b)  and  (9c)  will  next 
be  approximated  asymptotically  for  the  shadowed  portion  of  the  transition 
region  in  part  A of  this  section;  a corresponding  development  for  the  lit 
portion  of  the  transition  region  will  follow  in  part  B.  In  the  subsequent 
analysis,  <j>  will  be  restricted  to  0<<t><ir  for  convenience;  the  results  for 
— tt <4><0  may  be  readily  obtained  via  symmetry. 


Field  Analysis  for  the  Shadowed  Transition  Region 


(2) 

One  may  rewrite  (9b)  by  incorporating  (10)  and  by  replacing  HVx(kp) 
in  the  integral  over  cj  with  its  integral  representation. 


H^lko)  -1  Jde  e-3k»stn6tjv(,)B 


to  obtain 
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(12) 


I % ejkpcos|^  + 


-jkp[sinB-  ^ e] 


where  cr  is  shown  in  Figure  3.  The  integrand  in  (12)  possesses  a saddle 
point  at  e=05=cos"l  a/p,  and  a pole  at  B = 3D  = \p.  For  fixed  values  of 
p and  a,  8S  is  a constant  which  is  alternatively  given  by  3S  = <j>s  - -n/2, 
where  4>c  is  the  value  of  <p  at  the  shadow  boundary  (SB])  as  shown  in 
Figure  2.  The  angular  extent  of  the  shadow  zone  (for  0<<j><ir)  is 
thus,  Bp>6s  and  e>0  in  the  shadow  region,  where 

e = * - 0S  = ep  - es  = ♦ - ♦s.  (13) 


SD  P-J 

Im  J3  1 

kit 


Re/3 


f(M.I 

SHADOW  REGION 
FOR  0 <4><tt 


Figure  3.  Location  of  B§,  Bn  and  SDP  for  the  shadow 
region  analysis. 


The  shaded  areas  in  Figure  3 denote  the  regions  where  e-jkp^-sin^"  in  (12) 
becomes  vanishingly  small  for|lms|-*”  ; thus,  one  may  deform  the  path 
cr  in  (12)  into  the  steepest  descent  path  (SDP)  through  3S  as  indicated 
in  Figure  3.  Since  0§<3n  in  the  shadow,  the  pole  at  Bp  is  crossed  in 
this  contour  deformation;  therefore,  it's  residue  given  by  -eJ *cpcos<Ph( ) 
must  be  included  in  the  evaluation  of  the  integral.  H(A)  is  a 
step  function  which  is  unity  when  A>0,  and  zero  when  A<0.  Next,  the 
exponent  in  the  integrand  of  (12)  may  be  approximated  by  it's  three  term 
Taylor  expansion  about  S=Ss  for  large  kp  sin6s  to  obtain 
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•; 


* ejkpcos^  _ gjkpcos*  HU_^}  _ 


;-jk(s+ae) 

2ttj 


f 


de 

SDP 


where 


2 2 2 

s = p - a ; p sines  = s. 


J r 

e'eP 

(14) 


( 1 5a ; 1 5b) 


Note  that  the  first  two  terms  in  (14)  cancel  in  the  shadow  region.  This 
leaves  the  integral  over  the  SDP  which  may  be  directly  expressed  in 
terms  of  the  well  tabulated  Fresnel  integral  after  employing  the 
transformation 


-J 


v = e 


(16a) 


which  maps  the  SDP  onto  the  real  u-axis  (-«*y<").  Thus,  for  <#>><}> 


-jk(s+ae)  r Q-ksy 
‘l  * - — \d|J 


2 -jj-jkae 

F[kLa]  S— - , (16b) 


y-yp  e 


/s 


where  u is  the  value  of  u when  3*»e  , and 
P p 


F[kLa]  = 


= 2J|/kL*|  ejkLa  Jdx  e"JT  * 
|/kLa| 


(16c) 


A plot  of  F[kLa]  is  illustrated  in  [Kouyoumjian  and  Pathak,1974].  The 
physical  significance  of  the  geometric  quantities  0 and  s is  shown  in 
Figure  4.  Turning  next  to  the  evaluation  of  I2  in  (9c)  for  the  shadow 
region,  it  is  observed  that  the  integrals  in  (9c)  may  be  simplified  by 


employing  the  Debye  approximation  of  (11)  for  H^7(\(kp)  in  the  integrands. 
v(t)  = ka  + mx  * 0(ka)  in  the  transition  region, 'and  the  major  contri- 
bution to  the  integrals  is  for  t small;  hence,  one  may  use  sin  y % s/p  in 
the  Debye  approximation  to  obtain, 


H\/Mkp)  * 


-jks  + j J j (ka+mOs, 

s ' » 


Employing  (18)  and  (7c)  in  (9c)  leads  to 


where  P$(£)  are  Pekeris'  caret  functions  [Logan,  1959; Bowman  et.  al.,1969] 
defined  by 


PSU)  - 
h ^ 


= f StLsL  e'J{T  dT, 


_oo 


5^o(t) 


and  the  parameter  £ in  the  shadow  region  is 


£ = me  (£  > 0 in  shadow). 


The  subscripts  s and  h in  (20)  correspond  to  the  TMZ  and  the  TEZ 
cases,  respectively.  Plots  of  Ps(£)  are  illustrated  in  Figures  5a 

h 

and  5b  in  terms  of  related  functions  p*(£)  and  q*U).  Combining 
the  results  of  (19)  and  (16b)  according  to  (9a)  allows  one  to  write  u 
at  a point  P in  the  shadow  region  as 


j—  “J  ^ 

u(ps)  * “ — tr-  (l-F[kLa])  + P (?) 

2C  mt 
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Figure  4.  Ray  paths  for  the  shadow  region. 

The  above  result  may  be  interpreted  geometrically  as  a field  which 
propagates  to  P§  along  the  path  Qi Q2ps  a^ter  being  launched  by  the  in- 
cident wave  at  Q-j , as  shown  in  Figure  4.  Q2  is  the  point  of  tangential 
shedding  of  the  diffracted  field  from  the^surface.  In  tjje  deep  shadow 
region,  i.e.,  far  from  SB],  £»0  and  F[kLa]  1 since  kLa  becomes  large; 

hence,  for  5>>0  only  the  term  in  (22)  is  significant.  Furthermore,  for 
C»0,  the  integral  for  P^  in  (20)  may  be  replaced  by  a rapidly  con- 
vergent residue  series  [Logan, 1959] 


PSU) 


Jf«ne 

e _ j e _e 


57T 
"J  6“ 


/T  n 2[Ai'(-qn)]2 


e»o 


rJ? 


1 *= 


t 5tt 

i ir  — J 6 

J F 5qne 
e e 

T~ 


n 2qn[Ai(-qn)] 


for  TMZ  case 


(23) 


for  TEz  case 
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q*(*) 


The  Miller  type  Airy  function  is  given  by  Ai ( t)  = V ( and  A1'(t)  * 
d/dx  Ai(x).  The  parameters  qn  and  q-  are  defined  by  Ai(-qp)  = 0,  and 
Ai'(-qn)  = 0,  respectively  where  n = 1,2,3, •••;  their  numerical  values 
are  tabulated  in  [Logan, 19591.  Thus,  far  from  SB],  i.e.,  exterior  to  the 
shadowed  transition  region,  (22)  reduces  via  (23)  to  the  Keller  surface 
diffracted  ray  field  given  by 


Dp  and  afj  are  the  Keller  diffraction  and  attenuation  coefficients  for  the 
nth  surface  ray  mode  [Keller, 1956]  which  to  first  order  are  given  by 


s-,2 


[D“] 


m e 


[A1'(-qn)]2 


h 2 

[Dh„]  - 


m e 


■i  w 

-J  T7 


rzar  -= 


qnCA1(-qn)]2  ’ 


a 


S 

n 


a 


j 

m e 


IT 


The  superscripts  s and  h in  (24)  have  the  same  meaning  as  the  subscripts 
in  (20);  also  uMQi)  is  the  value  of  the  incident  field  at  Qi  which 
is  unity  in  this  case.  The  geometrical  interpretation  of  (24)  is  the 
same  as  that  in  Figure  4. 


CW  - 

Finally,  the  term  u corresponds  to  a residue  series  [Bowman  et.  al., 
1969]  which  represents  the  field  of  multiply  encircling  GTD  surface  ray 
modes  given  by 


ucw(P$)  -v  ui(Q1) 


s 

ih/ 


■(a>jk)(t+*T) 


1 1 TO  >•  TO) 

i=l  n 


+ ui(Qj ) 


,-jks 

7f“ 


■£,  

- _ k -(«5+jk)(f+[t-l]T)  J 
l l Dp(Qj)e  Dj(Q^) 

=1  n 


;-jks' 

Ss7 


, 0<  4><tt  ; 


— ' (25) 

where  T = 2ira;  t'  = arc  length  from  Qj  to  Q£;  and  s'  = distance  from  Q£ 
to  Ps.  The  first  term  in  (25)  is  the  field  of  surface  ray  modes 
launched  at  Q-j  which  diffracts  from  the  surface  at  Q2  to  arrive  at  Ps 
after  encircling  the  cylinder  t times  in  the  counter-clockwise  sense. 
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The  second  term  represents  the  field  <j)f  surface  ray  inodes  launched  at  Qj 
(see  Figure  4)  which  diffracts  from  Q2  to  Ps  after  encircling  the  cylinder 
U-l)  times  in  the  clockwise  sense.  For  sufficiently  large  ka,  only  n=l 
in  the  second  term  of  (25)  is  significant  and  the  remaining  contribution 


is  negligible. 
(22)  and  (25) 


The  total  field  at  Pc  in  the  shadow  zone  is  the  sum  of 


B.  Field  Analysis  for  the  Illuminated  Transition  Region 


If  one  directly  employs  the  result  given  in  the  previous  section 
for  the  shadow  region  to  calculate  the  field  u at  a point  P|_  close  to 
SB-j  in  the  lit  region,  then 


u(P|_) 


uMPl)  - in 


,-jkt 


(1 -F[kLa] ) + P^(sj 


,-jks 


t = ae  < 0. 


(26) 


Since  4><4>s  and  Pq<es  in  the  lit  region  H(<f>-<f>s)  = 0 in  (14).  Also, 

(19)  is  unchanged  except  that  5,t  and  0 now  become  negative  in  the  lit 
region.  The  second  term  in  (26)  may  be  interpreted  as  a scattered  field 
which  after  being  launched  at  Qi  propagates  backwards  from  Q]  to  Q3 
before  being  shed  tangentially  from  Q 3 to  P^  as  illustrated  in  Figure  6. 


Figure  6.  Reflected  and  pseudo  ray  paths  for  the  lit  region. 
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It  is  noted  from  Figure  6 that  this  "pseudo  ray"  propagation  path  is  dif- 
ferent from  the  geometric  optical  reflectea  ray  path  which  satisfies  the 
generalized  Fermat's  principle.  The  result  in  (26)  is  quite  accurate 
very  close  to  SB];  however,  just  as  the  result  in  [Fock , 1 951 ] and 
[Wait  and  Conda,1959],  this  result  also  does  not  reduce  to  the  geometric 
optical  ray  field  far  from  SB],  An  approximation  for  H(2J(kp)  which  is 
different  from  that  in  (18)  appears  to  be  necessary  to  adequately 
approximate  the  field  behavior  in  the  deep  lit  region.  The  following  ap- 
proximate procedure  is  adopted  to  achieve  this  goal.  First,  H(?^(kp) 

in  I2  of  (9c)  is  replaced  by  its  large  argument  approximation  (v(x)«kp) 
such  that 


l2 e 


-j(kp-  t ) 


4 1 

jm 


} ^o(t) 


c2  4 2 


$W-|(t)  jv(x)^  - 


f Qw ,(x 
+ y dx  —1 

i <Wt) 


c1  W"2V 


One  next  defines  the  following 


£ - * = 2y  ; 5'  - -2m  siny 


(28a; 28b) 


where 


(ar)  < K 


In  the  transition  region,  may  be  assumed  to  be  small  in  comparison  to 
m which  is  taken  to  be  large.  One  may  approximate  v(x)(n/2  -ip)  in  the 
exponent  of  (27)  via  (29)  as 

2" 

v(t)  (ip  -y'j  = (ka  + mx)2y  % -ka  jjj-  + ^ — x5'  . (30) 

m 
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Incorporating  this  result  into  (27)  yields 


. Mil 

5[w2(t) 


Referring  to  Figure  6,  it  is  seen  that  in  the  far  zone  (kp>>ka) 


$ % 201 ; p % 


i , in  amplitude 
a + a cose1,  in  phase 


2y  = it  - |<j>  | % n-  20 1 


(32a;32b) 


(32c) 


respectively.  Thus,  I2  may  be  expressed  in  terms  of  the  reflected  ray 
coordinates  as 

i i (e'13  r.i  w *"1 

rr  -jk(A+  a cose  ) j2ka  cose  -j  v J J 

I2'-Jp-e  me  ; 


5'  = -2  m cose1  (via  (28b)  and  (32c)).  (34) 

Turning  next  to  the  evaluation  of  I]  of  (9b),  one  begins  with  its  integral 
representation  in  (12) . The  contour  cb  is  deformed  into  the  SDP  through 
the  saddle  point  Bs  as  in  Figure  3;  the  pole  at  is  n°t  crossed  in  this 
contour  deformation  since  Sp^c  is  in  the  lit  region.  To  be  consistent 
with  (27),  I]  will  also  be  evaluated  in  the  far  zone;  thus,  es=cos_1a/pW2 
for  kp>>ka  so  that  a three  term  Taylor  expansion  about  6S  for  the  exponent 
in  (12)  yields 
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(35) 


I-j  % u1(PL)  - 


jka(f-  *) 


rJkp  U i 


kp[|<6-  J)  ] 


J 8-6 
SDP 


The  integral  in  (34)  may  also  be  directly  expressed  in  terms  of  the  F 
function  of  (16d)  which  contains  the  Fresnel  integral  as  done  previously 
for  the  SDP  integral  of  (14)  for  the  shadow  zone. 


-j(kp+  J)  jka  2y 

Ii  * uV.)  - £ — ; Z—z F[kp-2y2] 


Incorporating  (29),  and  the  far  zone  approximations  for  p and  <p  of  (32) 
into  (35)  yields  the  following  result. 


, e-JT  J2ka  cose'-j  -Jkt^cose1) 

I,  ^ u^Pi  ) + m e 2 F[kL'a']  — 

1 L vCTC  SI 

kp  » ka,  (37) 

2 

L'  = 1 ; a'  = = 2 cos2ei  (with  S'  as  in  (34)).  (38a;38b) 

2nr 

Combining  (33)  and  (37)  according  to  (9a)  gives  u(PL)  in  the  far  zone  as 


u(PL)  ~ u1  (PL)  + un(QR)  <- 


trn 


which  for  later  convenience  is  re-expressed  as 


(11)1  -J  5* 


(l-F[kL'a' ]) 


(39a) 


U(PL)  * ui(PL)  + ui(QR)  -JT 


+ Ps(6') 


jliilif  "J  t 

jr!  h — (l-F[kL'a']) 

4 2C'/7 


?y  -jk* 

- — P 

l 


(39b) 


where 


u1 (QR)  = incident  field  at  QR  (see  Figure  6)  = e^ka  cose  ; 
pY  = a cSse  = reflected  ray  caustic  distance. 


(39c) 


Far  frgm  SBi  in  the  lit  region,  £'<<0,  and  hence  F[ kL ' a ' ] 1;  further- 

more, PsU'j  may  be  asymptotically  approximated  as  [Logan, 1959; 

Logan  and  Yee,1962] 

j 

hu')',±v¥-  e ~yr  n+OU'"3)];  ?'«0.  (40) 


Thus,  incorporating  (40)  into  (39b)  for  the  deep  lit  region  U'«0)  yields 


u(PL)  * u1  (PL)  + ^(Qr)  Jp"e'jki  , 


)™zl 

for  the  \ > case.  (41 ) 


The  field  ubU  of  the  geometric  optical  incident  and  reflected  rays 
represents  an  accurate,  first  order  asymptotic  field  approximation  within 
the  deep  lit  region,  and  it  is  well  known  to  be 


uG0(Pl)  * u\pl)  + ui(QR)  Rs  J ^ — e' 


r 


{R  =-l  for  TM  

R -+1  for  TE  case 
h z 

it  is  apparent  Trom  tnat  tne  result  in  ^hij  is  tne  Tar  zone  geometric 

optics  field,  since  the  reflected  ray  divergence  factor j3'Y/(?'Y+t)  of  (42) 

reduces  to  4^/*  in  the  far  zone  (t»pY).  The  result  in  (39)  for  the  lit 
region  has  essentially  the  same  form  as  (22)  for  the  shadow;  this  is  a 
desirable  feature  since  the  results  in  (22)  and  (39)  are  in  terms  of 
tabulated  functions,  and  they  admit  a simple  interpretation  which  will  be 
discussed  in  Section  V.  The  far  zone  result  of  (41),  and  hence  (39)  may 
therefore  be  heuristically  generalized  to  the  near  zone  (which  is  exterior 
to  the  close  vicinity  of  the  surface)  according  to  ray-optics,  by  replacing 

the  far  zone  ray  divergence  factor  with  it's  near  zone  value 

where  «■  must  now  represent  the  near  zone  reflected  ray  distance  from 
Qr  to  P[_.  Thus,  the  near  zone  result  for  the  lit  region  is 


case. 


i(el)3J~JT 

u(PL)  - u1  (Pj_)  +u1(Qr)  - 


(l-F[kLl3f,])+Ps(cl ) 


. 2 e‘jkJl;  v = - 2m  cose1  < 0;  L'=*+£Y. 


It  is  observed  via  (22)  that  the  distance  parameter  L for  the  shadow  zone 
is  the  distance  s from  the  caustic  point  Q2  to  the  field  point  Ps.  The 
caustic  at  Q2  is  the  effective  origin  of  the  diffracted  ray  path  to  Ps. 

Thus,  keeping  this  interpretation  in  mind,  one  modifies  L'=*  in  (38a)  for  the 
far  zone  to  L ' = «.  +3^  in  (43)  for  the  near  zone,  since  A+f)*  is  now  the  distance 
to  Pi  from  the  virtual  caustic  of  the  reflected  ray  path  along  which  the 
field  u(P| ) is  assumed  to  propagate.  One  notes  that  L'  of  (43)  reduces  to 
that  of  (38a)  at  the  shadow  boundary,  and  also  in  the  far  zone 
(where  t»^).  The  value  of  L'  is  significant  only  near  the  shadow 
boundary  where  L'  of  (43)  and  (38a)  are  approximately  the  same;  hence  one 
could  use  L'  of  either  (43)  or  (38a)  without  noticeably  affecting  the 
numerical  calculations.  The  expression  in  (43)  now  properly  reduces  to 
(42)  in  the  deep  lit  region.  The  total  field  at  Pi  is  the  sum  of  (43) 
and  ucw  of  (25)  with  Ps  replaced  by  P[_;  however,  u^w  is  negligible  in 
the  lit  region  for  large  ka. 

It  will  be  shown  in  Section  V that  (43)  for  the  lit  region  and  (22)  for 
the  shadow  region  are  equal  at  the  shadow  boundary,  thereby  making  the 
asymptotic  solution  presented  in  this  paper  continuous  from  the  deep  lit 
region  to  the  deep  shadow. 
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III.  GENERALIZATION  TO  THE  SMOOTH  CONVEX 

CYLINDER  CASE 

The  uniform  results  for  the  circular  cylinder  in  (22)  and  (43)  may 
be  readily  generalized  to  treat  the  scattering  by  a convex  cylinder  of 
slowly  varying  curvature  in  the  usual  manner,  by  utilizing  the  local 
properties  of  propagation,  scattering  and  diffraction  of  waves  at  high 
frequencies.  Thus,  one  treats  each  point  on  the  convex  cylinder  as  if 
it  were  locally  on  a circular  cylinder  of  the  same  radius  of  curvature, 
i.e.,  the  radius  (a)  in  the  solution  of  the  previous  section  is  replaced  by 
pg  where  pg  is  now  the  local  radius  of  curvature  of  the  convex  cylinder. 
Specifically,  the  generalization  of  (22)  and  (43)  to  the  case  of  plane 
wave  scattering  by  a smooth,  perfectly-conducting  convex  cylinder  of 
variable  curvature  is  indicated  below  for  the  shadow  region  in  part  A, 
and  for  the  lit  region  in  Part  B. 

A.  Shadow  Region 

The  field  u at  Ps  in  the  shadow  zone  of  the  convex  cylinder  becomes 


u(P$  * u1^) 


-/m(Q1)/m(Q2) 


e'Jkf  o-FtBflHyo: 


,-jks 


(44a) 


5>0 


with 


r . r 

e = \ df  ; t-j  df  ; m(f)  = 

Q1 


kp  „(f ) 


1/3 


(44b,44c,44d) 


and 


L = s 


a = £2/[2  m(Q1)m(Q2)] 


(44e;44f ) 


The  points  Qi  and  Q2  are  indicated  in  Figure  1;  t'  denotes  any  point 
between  Q-j  and  Q2  on  the  cylinder  and  s is  tl]e  distance  from  Q2  to  Ps. 

It  is  noted  that  m in  (22)  is  replaced  by/m^Qi ) /mTQ2 ) in  (44a)  to  be 
consistent  with  reciprocity,  since  the  point  of  grazing  incidence  Qi  and 
the  point  of  diffraction  Q2  on  the  convex  surface  in  general  possess 
different  radii  of  curvatures  Pg(Q-| ) and  Pg(Q2).  respectively.  This  sym- 
metrical splitting  of  m is  also  essential  for  preserving  the  uniform 
property  of  the  result  so  that  deep  in  the  shadow  region  where  ?»0,  (44a) 

reduces  to  the  GTD  result  of  (24)  (except  that  (a[^  t)  therein  is  replaced 


by 


f2  s 

Jn 


(t')dt'  for  the  convex  cylinder). 
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B.  Lit  Region 

The  field  u at  PL  in  the  lit  zone  of  the  convex  cylinder  becomes 


u(PLhui(PL)+ui(QR) 


-i  111 ) ( ~Jt 

12  (l-F[kL'a'])+P  (5 

i26'/7  i 


■l 


g-e’^1;  5'  < 0, 


(45a) 


with 


V = - 2 m(QR)cosei  ; m(QR)  = 


W 


il/3 


and 


L'  = l + p*;  a'  - 2cos2ei. 


pq(QR)cosel 

; pY  = -2—^ , 

(45b;45c;45d) 


(45e;45f) 


The  point  of  reflection  Qr  is  shown  in  Figure  1;  also  £ is  the  distance 
from  Qr  to  P|_  as  before.  Deep  in  the  lit  zone  where  ,T<<0»  (45a)  reduces 
to  the  geometrical  optics  result  of  (42)  except  that  p is  now  given  more 
generally  by  (45d)  for  the  convex  cylinder. 


IV.  FIELD  IN  THE  CLOSE  NEIGHBORHOOD  OF  A SMOOTH  CONVEX  CYLINDER 


While  the  uniform  result  in  (44)  and  (45)  for  the  scattering  of  a 
plane  wave  by  a smooth,  perfectly  conducting  convex  cylinder  is  valid  in 
the  far  and  near  zone  of  the  surface,  it  is  not  valid  in  the  close 
neighborhood  of  the  surface.  The  region  of  validity  of  (44)  and  (45) 
excludes  the  boundary  layer  and  it's  vicinity  in  the  shadow  region,  and 
it  excludes  approximately  the  same  distance  from  the  surface  in  the 
lit  region.  The  field  in  the  close  neighborhood  of  the  surface  is 
available  elsewhere  in  terms  of  the  canonical  Fock  integral  [Logan  and 
Yee, 1962; Ivanov, 1971 ] which  is  very  complicated  for  numerical  calculations. 
In  this  section  a Taylor  series  approximation  for  the  canonical  Fock 
integral  is  obtained  to  describe  the  field  in  the  close  neighborhood  of 
the  surface.  This  approximation  is  very  convenient  for  numerical 
calculations. 

When  the  field  point  is  extremely  close  to  the  surface  of  a circular 
cylinder  such  that  k(p-a)<<ka,  one  may  then  approximate  J^^j(kp)  and 

Hi(l)(kp)  in  t8®)  and  (8b)  by 


Jv(kp) 


0) 

i(2), 


(m/IT)  ’v(x-h);  H^(kp)  % + j(m/rT)  ^w-|(x-h) 

2 


(46a); (46b) 


where 

h = nf^kd;  d = p - a.  (46c);(46d) 

Incorporating  the  above  results  into  (8a)  and  (8b)  essentially  yields 
an  expression  for  the  field  u in  terms  of  Fock's  canonical  integral  as: 


u % - 


-jkaifi 


( dT[V(x-h)  - w,(x-h)]e"^ZT;  Z = nty  (47a) 

y . four)  2 


or 


r""je  $w  (t) 

u*li_e-jkaM  dx[w,(x-h)  - 5-1 w?(x-h)]e-jZT;  Z=m*.  (47b) 

2^  J_je  1 ?fW2(x)  2 


with  tp  = |<j>|  - tt/2  as  before. 
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The  above  Integrals  for  the  circular  cylinder  case  will  be  approximated 
in  the  close  vicinity  of  the  surface  by  a Taylor  series  in  direct  powers 
of  the  normalized  distance  h;  these  results  will  then  be  generalized 
to  the  convex  cylinder  case  as  in  section  III. 


The  ipO  case  will  be  considered  first,  to  be  followed  by  a similar 
treatment  for  the  ip<0  case.  Let  0<<f><ir  as  in  section  II;  the  results  for 
tt< <#><2tt  may  be  directly  obtained  via  symmetry.  For  it  is  convenient 
to  employ  the  representation  in  (47b).  Thus,  for  h small,  one  may  repre- 
sent w,(x-h)  in  (47b)  by  a Taylor  series. 


w^(x-h)  ■ w^(t)  - w-j(x)h  + tw^(t) 
2 2 2 2 

+ [2w^(t)  + t^W^(t)]  jq  - 


[w-j(t)  + ^(t)]  £- 

2 2 

5 

[4tW-|(t)  + T^(x)]  yrjjj-  + 0(h®) 


(48) 


where  the  primes  on  w-j(t)  denote  differentiation  with  respect  to  t.  In- 

2 

corporating  (48)  into  (47b)  allows  one  to  express  u in  terms  of  Fock 
current  functions  and  it's  derivatives. 

u(p)|ipo  * e_jka,,J^g(z)  + lr  3'(z>-  g(z)_  {£g"(z7jH)(h6); 

for  the  soft  or  TMz  case  ($=1  case) 

(49a) 

u(p)|^>o  ^ e"jka^(z)  +jr9'(z)  -3T9(z)  - JfT  9"(z)  -$T“  g'(Zj] 

+ O(h^);  for  the  hard  or  TE^  case  (?[  = 3/9x  case).  (49b) 
g(D)  and  g(D)  are  the  Fock  functions  [Fock, 1946;Logan, 1959]  defined  by 


-je 

dx 

-Je 


e-jDt 

wjrn’  • 


(50a;50b) 


The  values  of  q(D)  and  g(D)  are  extensively  tabulated  [Logan, 1959], 

The  primes  on  g(D)  and  g(D)  in  (49a;b)  denote  differentiation  with  respect 
to  the  argument  (D).  The  derivatives  of  % and  g may  be  evaluated  numerically 
from  the  interpolated  values  of  g and  g,  respectively.  In  obtaining 
(49a;b),  use  has  been  made  of  the  well-known  relationships  w^'(x)  = tw^(x); 

2 2 
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w](t)w2(t)-w1(t)w^(t)=-2j;  and 


p e-JOT  Tnf(t)  . j"  il  p e-JD,  t(T), 


-“-Je 


where  f(x)  represents  some  combination  of  the  Airy  functions.  It  is 

noted  that  g(Z)e"^ka,Jj  and  g(Z)e"3ka<^  are  directly  proportional  to  the 
electric  current  density  induced  on  the  cylinder  surface  by  the  in- 
cident plane  wave  [Fock,1946].  In  the  present  case,  these  currents  exist 
at  Pn  which  is  the  normal  projection  on  the  surface  of  the  field  point  P 
exterior  to  the  surface.  For  <p0,  the  point  Pr  lies  in  the  shadow  zone. 
The  distance  from  Pw  to  P is  d(=p-a);  furthermore,  the  unit  outward 

surface  normal  at  Pfj  is  p = n = PfjP/  P^T*  . When  the  point  P^  on 
the  surface  lies  in  the  deep  shadow  (^>>0  and  hence  Z»0),  the  Fock 

currents  g e-jka^  and  g reduce  uniformly  to  the  GTD  current 

representation  [Pathak  and  Kouyoumjian,1974]. 

For  ip<0 , it  is  convenient  to  begin  by  re-writing  (47a)  in  terms  of 
a new  parameter  V as  employed  for  the  Fock  currents  by  Logan  (1959). 


m simj>  = -m  cos <(>;  \p 


(51a, 51b) 

The  approximation  in  (51b)  is  valid  within  the  shadow  boundary  transition 
region  for  large  m.  Thus,  (47a)  becomes 


u(p)|.<0  * eJka  cos*  “ 


j(Z')3/3 


-•-je 


(dx  V(x-h)  - (T_h)  e"jZ'T; 

l 5fwo(x)  2 


(52a) 


u(P)  .<fl  % e3ka  cos*  e-jhZ’ 


■j  (z1  r/3 


S*  St, 


-•-je 


(52b) 


In  obtaining  (52a),  one  employs  e"^T  % e”3^  T and  e”3^  % e^a  cos<#>~j ) /3. 
Furthermore  employing 


”-Je 

k id;' 


V(T-h)e 


-JZ't  _ -jhZ1  J( Z')3/3 


= e 


for  ii>< 0 in  (52a)  leads  to  (52b).  Next,  the  v^r-h)  in  (52b)  may  be  expanded 
for  small  h as  in  (48)  to  arrive  at  an  expression  similar  to  that  in  (49a;b) 
for  u in  terms  of  the  Fock  currents  and  it’s  derivatives. 


u(l>)|,<0  % ejka  “=♦ 


+ e-j(Z’)3/3jhg(Z,)t  |£g'(Z')  -^g(Z')  -l£g"(Z')j 
+ O(h^),  for  the  soft  or  TM  case  (fy=l  case). 


(53a) 


u(p)|*<0*  eJka  C0S<P  e"jhZ'  ( J hZ ' )nj 

+ e-j(z,)3/3|g(Z')  + *jy-  g ' (z ' ) - ^g(Z')  - j£g"(Z')  - 
l'(Z')j 


j 5T  g' 


+ 0(h®),  for  the  hard  or  TEZ  case  (ft  = |^-  case) 

(53b) 


For  ip<0 , the  point  on  the  surface  is  in  the  lit  zone,  where  is  the  normal 
projection  on  the  surface  of  the  field  point  P,  as  indicated  previously. 

The  Fock  currents  eJka  “s-MU'W  g(Z')  end  eJka  ™s*-j(Z'P/3  g(z, X.jMP|))) 
at  Pfj  which  are  induced  by  the  incident  plane  wave  reduce  uniformly  to  the 
geometrical  optics  currents  when  Pm  moves  into  the  deep  illuminated  region 
(Z'«0). 
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The  results  in  (49a;b)  and  (53a ; b)  may  now  be  generalized  as  In 
section  III  to  obtain  the  field  u in  the  close  neighbourhood  of  a smooth, 
perfectly  conducting  convex  cylinder  illuminated  by  a plane  wave.  Thus, 


u(P)  % T.|+T2\hg(D)+  g'(D)  - ^y—  g(D)  - {jy  g"(D)>  + O(h^), 


for  the  soft  or  TMZ  case. 


(54a) 


u(P)  % U0e‘j°  T,  + ^/g(D)+  ^g'(D)  - g(D)  - gH(D)-  ^ g'  (D)j  +0(h6) 

— rr.L\ 


for  the  hard  or  TEZ  case. 


where 


(Q-j } 


k \ dt,! 


;n 


0,  9 


-mCP^Jcose1 


for 


P„  in  shadow 
" zone 


P.,  in  lit 
N zone 


(55a;55b;55c) 


e.jhD  . J Mi!  (jhD)" 


fpg(PN^ 

T J W 

2 i — 


i D 

"J  T 


; for 


P.,  in  shadow 
N zone 


Pw  in  lit 
. N zone 


(55d;55e) 


The  arc  length  integrals  in  (55b)  and  (55c)  are  evaluated  along. the 
surface  geodesic  from  Q]  to  Pfj  in  the  shadow  region;  whereas,  is  the 
acute  angle  between  the  incident  ray  and  the  surface  normal  to  Pr  in  the 
lit  region.  While  the  definitions  for  u0,  o and  D for  the  convex  cylinder 
follow  fairly  directly  from  the  generalizations  of  the  quantities  in  (49a, b) 
and  (53a;b)  for  the  circular  case,  the  introduction  of  [pg(PN)/Pg(Ql )]^6 
in  To  for  Pfj  in  the  shadow  is  required  to  preserve  the  uniform  nature  of 
the  Fock  currents  in  the  shadow  zone  of  the  convex  cylinder  [Pathak  and 
Kouyoumjian,1974]. 
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V.  NUMERICAL  RESULTS  AND  DISCUSSION 


A uniform  asymptotic  result  for  the  near  zone  field  exterior  to 
a perfectly  conducting  convex  cylinder  which  is  illuminated  by  a plane 
wave  as  shown  in  Figure  1,  is  given  in  (44a)  and  (45a)  for  the  shadow 
and  lit  zones,  respectively.  These  results  may  be  employed  via  reciprocity 
to  directly  obtain  the  radiation  patterns  of  electric  or  magnetic  line 
sources  in  the  near  zone  of  the  cylinder.  A few  numerical  results  for  the 
radiation  patterns  of  line  sources  near  circular  and  elliptic  cylinders 
which  are  based  on  (44a)  and  (45a)  are  presented  in  Figures  7b,  7c  and  7d. 

The  accuracy  of  this  uniform  result  is  confirmed  by  the  excellent  agreement 
obtained  between  the  present  solution  and  the  eigenfunction  solution  for 
the  circular  cylinder  case.  Similar  agreement  is  observed  for  the  elliptic 
cylinder  examples  in  which  case  the  present  result  is  compared  with  a moment 
method  solution.  The  results  in  (44a)  and  (45a)  are  not  valid  in  the  close 
vicinity  of  the  surface;  hence,  the  results  in  (53a)  and  (53b)  may  be  used  to 
calculate  the  field  in  the  close  neighborhood  of  the  surface.  In  particular, 
the  latter  results  are  used  to  calculate  the  far  zone  fields  of  line  sources 
in  the  close  vicinity  of  a circular  cylinder  via  reciprocity  as  shown  in 
Figures  8a  and  8b.  It  is  seen  that  these  results  tend  to  blend  with 
those  obtained  via  (44a)  and  (45a)  thereby  providing  an  overlap  region  for 
these  two  representations.  It  is  noted  from  the  above  pattern  calculations 
that  the  uniform  result  of  (44a)  and  (45a)  makes  the  leading  term  of  the 
total  high  frequency  field  continuous  at  the  shadow  boundaries.  On  the 
other  hand,  the  use  of  the  usual  ray  optical  or  6TD  solution  would  have 
produced  a discontinuity  in  the  field  pattern  at  the  shadow  boundaries. 

One  may  verify  the  continuity  property  of  the  uniform  result  as  follows. 

Near  a shadow  boundary,  ? and  V must  both  approach  zero  since  the  arc 
length  t 0,  and  e1  ir/2.  Furthermore, 


F(o) 


/iro"  - 2ae 


jfl  3(f  + 

e 


a) 


as  o 0, 


(56a) 


and  one  may  re-express  Pc  as 


Ps(o)  = 


s 

p*(°)l  ! 

[q*(a)(  2/JTa 


■i" 

■J4 
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(56b) 


VIA  EQS  (44a)  ft  (45a) 


figure  7c.  Total  field  surrounding  the  cylinder  of  Figure 
7a  excited  by  a magnetic  line  source. 


Incorporating  (56a)  and  (56b)  into  (44a)  and  (45a)  yiVds  exactly  the  same 
limit  value  as  the  field  point  approaches  the  shadow  boundary  from  either 
the  shadow  or  the  lit  side;  this  limit  value  is 


'V, 


1 

2 


u^P)  - m(Q]) 


TEz  case 


TMz  case 


where  *Sn  is  the  distance  from  Q]  to  P on  the  shadow  boundary  (SB).  Thus, 
(44a)  ana  (45a?)  yield  a continuous  total  field  at  the  shadow  boundary.  It 
can  be  shown  for  the  circular  cylindrical  case  that  the  total  high 
frequency  field  based  on  (44a)  and  (45a)  possesses  derivatives  along  the 
normal  direction  to  the  shadow  boundary  which  are  also  continuous  to 
0(1/asb)*  In  this  connection,  it  has  therefore  been  observed  that 
the  pattern  in  addition  to  being  continuous  is  also  relatively  smooth. 
These  observations  hold  also  for  the  smooth  convex  cylinder  case  provided 
the  local  radius  of  curvature  of  the  convex  cylinder  is  slowly  varying 
as  assumed  previously.  The  result  of  (53a)  and  (53b)  is  continuous 
across  the  shadow  boundary  as  may  be  easily  verified  by  inspection. 


While  the  uniform  result  of  (44a)  and  (45a)  pertains  to  a perfectly 
conducting  smooth,  convex  cylinder,  a similar  result  can  be  obtained  for 
a cylinder  with  an  impedance  boundary  condition.  However,  the  latter 
case  has  not  been  included  here  since  the  functions  which  are  now 

dependent  on  the  non-zero  value  of  the  impedance,  are  not  extensively 
tabulated. 


Further  extensions  of  the  uniform  result  of  (44a)  and  (45a)  to  treat 
the  general  three  dimensional  problem  of  the  scattering  of  an  arbitrary 
electromagnetic  ray  field  by  an  arbitrary,  smooth,  perfectly  conducting, 
convex  surface  will  be  presented  in  a future  paper  [Pathak,  et.  al . ] . 
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